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Abstract: We argue that heavy- ion collisions provide the best testing ground for mini- 
black hole physics as Mp ~ AGeV for the gravity dual of YM and give concrete evidence 
for a new extra dimension, that is visible only to the strong interactions. We analyse 
the process of production evolution and decay of the mini-black-holes by using recent 
results on gravity duals of YM. There are several novelties compared with the traditional 
story of black hole evaporation, including Bjorken scaling instead of sphericity, evaporation 
via bubble nucleation instead of the Hawking mechanism and lepton-poor final states. 
Multiplicities are estimated using shock- wave scattering techniques. It is argued that high- 
multiplicity/high energy pp collisions will also show similar characteristics of mini-black- 
hole production and decay. 
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1. Introduction 

The subject of the EPS talk was a review of holographic model of four dimensional YMs, 
Q , its implications for the finite temperature structure of the theory, ^ , and its predictions 
for the calculation of transport coefficients Q and the Langevin diffusion of heavy quarks 
in the quark-gluon plasma, A review of the model can be found in ||. 

What we would like to do here, is to explore the physics of the heavy-ion collisions in 
the dual gravitational description and argue that this process can be described in the dual 
gravitational description as the collision of two energy lumps that form a trapped surface 
that eventually leads to a horizon formation, and ensuing thermalization. The mini-black 
hole thus produced is unstable and it evolves classically until it decays via evaporation to 
a large multiplicity of final states. 

Such a picture was explored in the context of an AdS theory with a hard IR cutoff by a 
number of authors, || 0]. Since then, the structure of the bulk theory is better understood, 
and we will analyze the process in view of the new developments. 

Black holes have dominated the mysteries associated with gravity over the past decades. 
It therefore seemed like a golden opportunity when a production of mini-black-holes seemed 
credible, g, g at energies not far from the TeV region, if gravity had a lower characteristic 
scale. This provided hope that the mysteries of black-holes, namely their thermodynamics 
behavior, their entropy of extreme size, the Hawking radiation and the ensuing informa- 
tion paradox, could be studied at colliders. Experimenters could provide an experimental 
scrutiny that looks more promising that expecting indirect signals from astrophysical black 
holes that are "stars" these days of galaxy-core gravitational phenomena. 

Theories of large extra dimensions are still candidates for reality, however current ex- 
perimental constrains make the size of a higher-dimensional Planck scale recede further 
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from the currently accessible energy range, making the possibility of mini-black hole pro- 
duction look rather remote. 

Despite this unfortunate state of affairs, the gauge theory/string theory correspondence 
provided an alternative arena for mini-black hole production. QCD, the gauge theory 
of strong interactions is expected at large iV c to be described by a string theory. This 
theory seems elusive, as the string is very soft at high energies, but the expectation of a 
gravitational description at low energies is well founded, and a model that has the right 



properties and matches many YM observables has been proposed in Q and [10|. 

We will argue further on, that the heavy-ion collisions at CERN and LHC as well as 
the proton-proton collisions at LHC with high-multiplicity final states can be described by 
the formation and evaporation of a mini-black-hole, possible as the 5-dimensional Planck 



scale has been determined to be Mp ~ 4 GeV, [11]. 

Moreover, the gravitational theory has a single extra dimension, |j|], that is indirectly 
visible via the gravitational physics. This dimension may be expected to be a bit "fuzzy" 
as the the number of colors of the strong interactions N c = 3 seems far from the large- N c 
limit. However, as lattice simulations suggest (see figure [j] and [12]), SU(3) YM seems 



very close to SU(oo) YM. It is therefore fair to say that most of the quantum gravitational 
effects, controlled by are small. 

The process of black hole formation and evaporation in the dual theory of QCD has 
similarities with expectations for standard asymptotically flat black-hole production in 
gravity, like 

1. large total cross section, 

2. very large multiplicity events, but with few hard final particles 

3. suppression of hard perturbative scattering processes, 

4. Democratic primary decay (but only in the strong interactions). 

However, the process of black hole formation and evaporation in the dual theory of 
QCD differs in several respects from the processes that many authors advocated for LHC, 
namely 

1. The gravity in question is not the standard gravity, but the one of QCD, where the 
gravity is mediated by the 2 ++ glueball. Other fields, and most importantly the 
"dilaton" (0 ++ glueball) play also an important role. 

2. Instead of high sphericity events, there is approximate Bjorken boost invariance in 
the metric. 

3. They are BHs in asymptotically AdS space (almost) instead of asymptotically flat 
space. 

4. Their primary decay is not via Hawking radiation but via bubble nucleation. 

5. The ratio of hadrons to leptons is small. 
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The main reason for the differences is that the BHs thus generated are live in a curved 
space, and that the theory contains also other fields (in particular the dilatonic scalar). 



2. The static, translationally-invariant black holes of IHQCD 

These are the translationally invariant saddle-point solutions, of IHQCD, defined by the 
Einstein-dilaton Lagrangian, j|], 

S 5 = -MlNl [ d 5 x^ R-\{d4>) 2 + V{<t)) +2M„ 3 iV c 2 / d 4 xVh K. (2.1) 

J L ^ J JdM 

Here, M„ is the five-dimensional Planck scale and iV c is the number of colors. The last 
term is the Gibbons-Hawking term, with K being the extrinsic curvature of the boundary. 
The effective five-dimensional Newton constant is G5 = l/(167rMp N 2 ), and it is small in 
the large- A^ limit. 

Of the 5D coordinates {xi, r}j = o...3, Xi are identified with the 4D space-time coordi- 
nates, whereas the radial coordinate r roughly corresponds to the 4D RG scale. We identify 
A = with the running 't Hooft coupling N c gy M , up to an a priori unknown multiplicative 
factor. 

The dynamics is encoded in the dilaton potential, V(A). The small- A and large-A 
asymptotics of V(A) determine the solution in the UV and the IR of the geometry respec- 
tively. 

1. For small A, V(A) is required to have a power-law expansion of the form: 

V(X) ~ + + + ...), A^O. (2.2) 

2. For large A, confinement and linear Regge trajectories (m 2 ~ n) require: 

V(X) ~ At y/\ogX A -> 00, (2.3) 

The theory, upon tuning two phenomenological parameters describes well both zero- 



temperature spectra of glueballs as well as the thermodynamics at finite temperature, [11]. 
In particular at T = the ground state is described by a unique saddle-point solution. For 
< T < 

Tmini this is the only saddle point available to the system, upon compactification 
of the time-circle. We will call this solution the "thermal vacuum solution" as it describes 
the confining vacuum of large- A^ YM. For T > T m i n , there are generically two saddle point 
solutions, (beyond the thermal vacuum solution), that are black holes, Q as can be seen in 
figure ||. The one with the larger horizon, is a black hole with positive specific heat while 
the smaller one has a negative specific heat. These two black holes merge at T = T m i n . At 
T = T c > T m i n there is a first order phase transition to the deconfinement phase, described 
by the large black hole, [|j. At T > T c the theory is in the deconfined phase. A calculation 
of the thermodynamic functions in this phase, and its comparison to a recent high statistic 
lattice calculation is shown in figure |]. 

The basic characteristics of the translationary invariant black hole solutions of IHQCD 
are as follows 
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Figure 1: Left: normalized entropy as function of temperature for YM with various N c — 
3,4,5,6,8. Right: the normalized trace of the stress tensor (interaction measure). The yellow 
lines are based on the calculations of IHQCD. Data by M. Panero from [fl2|. 



The large black holes have a E(T), S(T) that matches the one in large N c YM. In 
particular E ~ T 4 , S ~ T 3 as T — > oo. The horizon position is rj^ rge 
large temperatures, where r 
heat is positive. 



y db\j 

is the position of the AdS boundary. The specific 



The small black holes are unstable, with negative specific heat. They are however 
nowhere near Schwartzschild black holes. Although like Schwartzschild black holes 
they are unstable, they have scalar hair that makes them have different properties 
as obvious from their thermodynamic functions. As T — > oo their horizon shrinks to 
zero size as r? ma11 ~ -r^- — . Note that r — > oo is the IR region of the bulk geometry. 

A QCD 

We also have 



S ~ V3 exp 



T 2 



A 2 

QCD 



E ~ V 3 Mp T exp 



A 2 
QCD 



OO 



(2.4) 



It is interesting that at large T, the small and large black holes satisfy the duality 
relation 

(A QCD rr^)(A QCD gc ) ~ 1 (2.5) 



which should be compared with the analogous relation in global AdS, (£ r h 



small 



)(£r 



large \ 



1. 



Using as parameter r^, we can put both the large and small branches on the same 
diagram. For < r m i n we are in the large bh branch while for > r m i n we are in 
the small bh branch. The two branches merge at r m i n at a common temperature Train- 
The functions E(r n ), S(r n ) are monotonic (decreasing) functions of r n . As functions of 
temperature, in each branch, E(T), S(T) are monotonically increasing functions. They 
also satisfy, Q 

Elarge(Ti) > E sma .\\(Tj) , Slarge(7i) > S Bma n(Tj) (2.6) 

for any Ti,Tj. 
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Figure 2: Typical plot of the black- hole temperature as a function of the horizon position rh, in 
a confining background. The temperature features a minimum T m i n at r m i n , that separates the 
large black-hole from the small black-hole branches. The boundary is at r = and the IR endpoint 

at r = oo 

3. Black holes in heavy-ion collisions 

Although the above properties are for translationary-invariant black holes, we can envisage 
finite configurations with such profiles. They will not be stable but will have thermody- 
namic characteristics that are similar. 

A particularly important question is what kind of black hole a heavy-ion collision 
creates. In a sense this is a question for the microcanonical ensemble, as there is a finite 
available energy E{ to create the black hole. The questions is, given the energy, what 
maximizes the entropy, as this will be the most likely outcome of the process. 

Noting that the energy density is a monotonic function of , there is a limiting energy 
density, namely e m i n = e(r m j n ) that separates large from small black holes. There is also 
another benchmark, namely e c = e(r c ) > e m i n at the critical radius, where the first order 
phase transition takes place. 1 

Therefore, if the energy density e« of the heavy-ion collision is bigger than e c we expect 
the formation of a large finite size bh. If e« < e m i n it is plausible that a small black hole is 
formed, although this is issue is rather complicated even in the case of AdS, |/], |l3| ]. The 
case e c > 6{ > e m ; n is murkier as in this case there are two possibilities. Either the creation 
of a large black hole in a subdominant saddle point, or the evolution of energy without 
thermalization, in the thermal vacuum solution. At N c — > oo the difference of the saddle 
point energies scales as N% and therefore no black hole is expected to form. However for 
SU(3) YM we might expect a potentially significant chance of thermalization. In view of 
the crossover behavior of real QCD, this regime is murky. However at RHIC and LHC, 
€i > e c . 

Given an initial energy density ej > e c , distributed in three volume V3 we expect 
that very fast it will form a trapped surface, [14, 18 1, a signal of thermalization and of 



developing a horizon. For V3 >> AqJ^ d as expected in heavy-ion collisions, a large black 

In large- N c YM the transition is strongly first order. SU(3) YM still has a (weaker) first order phase 
transition. Once massive quarks are included the transition becomes a crossover, but for our purposes the 
crossover is fast enough so that we t treat it a first order phase transition at the level of accuracy we are 
aiming at. 
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hole is formed with finite spatial extend roughly given by V3. Outside this extend, the 
saddle point solution is that of the vacuum solution. 

Because of the finite spatial extend this solution is unstable and begins to evolve with 
time. The main form of evolution is that it "falls" in the bulk geometry. On one hand 
there is a "surface tension at the spatial boundaries where there are big gradients, and the 
black hole expands in them. At the same time the radial horizon position recedes to the IR 
as expected from a "falling" energy distribution in the bulk. This keeps that horizon area 
approximately constant, and is therefore an isentropic evolution. The way this happens is 
that the 4d-profile remains invariant in comoving coordinates. It seems shrinking because 
it moves to larger values of r, where the metric is different. Therefore, because of the 
warping, the size in transverse space expands in physical units while the size in the radial 
direction contracts (in physical units). This is an adiabatic expansion in 4d QCD terms. 

In realistic heavy-ion collisions, there is deposition of about 20 TeV of energy (RHIC) in 
the midrapidity range during the first 0.1-0.5 fm/c with initial density about 15 GeV/fm 3 . 
From 0.5 - 1 fm/c there is ultrafast thermalization and entropy production. The free-fall 
phase is from about 1-6 fm/c during which the entropy is constant. Finally in the last 6-10 
fm/c this is an isothermal phase dominated by the transition. 

During the free-fall phase the Hawking temperature of the black hole is a decreasing 
function of time. It starts from e, = e(Tj) and it changes while keeping the entropy constant. 
As the entropy of the large black holes is that of YM, the gravitational and plasma pictures 
are identical. 

We may use the first law to calculate e'(T) = Ts'(T) and s' — = ^ with A the 
conformal anomaly, A(T) = e{T) — 3P. This relation can be integrated to 

s{T) = Sc T ,+T ^ A («)- , e(T) = + — j^A („)- + — (3.1) 

where T c is the temperature of the lrst order phase transition and S c is the transition 
entropy in the deconfined phase. We also obtain e c = e(T c ) = \T c s s + -f Below the 
transition we can approximate e ~ s ~ A ~ 0. Continuity of the free energy implies 
A c = T c s c , e c = A c . In Pure YM, 

T c ~ 240 MeV , e c = 0.31 T c 4 ~ (310) 4 MeV 4 , s c ~ ^ = (338) 3 MeV 3 (3.2) 

At the begining of thermalization the energy density is £j ~ 5.4 GeV//m 3 . (l/fm=200 
MeV). The system reaches the cross-over temperature T c after about t a d — 6 fm/c. 

For an adiabatic evolution dS = 0. if we define the linear dimension L as V3 = L 3 , 
and as the system expands with the speed of light we have ^ = Lo+ct • F rom figure |l] we 

find that s ^ c ) ~ 14 from where we find that Yf- = 14 or ~ 2.4. From this we find that 

s{T c ) Vo L Q 

ct ad ~ 1.4L , or L ~ 4.3 fm. 

When the black-hole temperature reaches T c , the solution becomes unstable. Unlike 
the standard Hawking picture of evaporation, the decay here proceeds by bubble nucleation 
where the confining vacuum is spontaneously appearing in several places seeded by inho- 
mogeneities in the distributions. This process is nearly isothermal, and can be in principle 
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described classically upon bubble nucleation. The decay time can be in principle estimated 
from the thermodynamic functions of the black hole. 

Small black holes on the other hand Hawking evaporate, in the standard fashion, but 
this process is not relevant for heavy-ion collisions. 

It is important to follow this black-hole evolution by solving the Einstein-dilaton equa- 
tions directly along the lines of [15, [R], 17, 18 1. 



4. Entropy production and particle multiplicities 

The collision of two energy distributions, as realized in heavy-ion collisions, can be ap- 
proximated at RHIC and LHC energies as the gravitational collision of two shock waves. 
Following Penrose and Eradley and Giddings we will search for a trapped surface that 
forms in the part of the space-time that is described by the superposition of two shock 
waves, before they start interacting. This trapped surface will eventually evolve to become 
a horizon. The area of this trapped surface therefore will give us a lower bound on the 
generated entropy during the collision, which will be converted to a multiplicity for the 
hadronic final state. Such an approach was developed recently for the AdS geometry in 



The metric and dilaton for such shocks can be found by solving the equations stemming 
from (|2.1|) with shock wave ansatz 



ds 2 = b{rf [dr 2 + dx i dx i - 2dx + dx- + ^(r^ 1 ,x 2 )5(x + )(dx + ) 2 ] , = (j)(r,x + ) (4.1) 
Compatibility of these equations implies that d+cj) = 0. The equations for $ are 

vi + 3^-d r + af] $ = -24 v 2 ± = didi , 4 = 8ttg 5 (4.2) 



b 

where we have introduced a stress-tensor source J++. 

The trapped surface and its boundary C, is determined by a function $ satisfying the 
differential equation and boundary conditions 

(a AdS . i -A(b(r)))(^-^) = , *|c = , Yl ( V ^) 2 !c = 8b 2 (4.3) 

j=l,2,r 

We may then compute the area that will provide a lower bound on the entropy 

1 f I 1 f r °i 

S > Sapped = l^T \J det \9AdS 3 \ dz ^ X S- = J Hrfx^r) (4.4) 

where the (generalized) curve C defines the boundary of the trapped surface S and where 
we have included the two sections of the surface associated with the two shocks. 

4.1 Shocks without transverse dependence 

This is the simplest case of a shock which corresponds to an infinite homogeneous nucleus 
on the transverse plane. Although the transverse distribution is unrealistic, such a simpli- 
fication allows for tractable calculations while it describes qualitatively various processes of 
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QCD [^7], 28]. As a first step, we begin our investigation using these kind of shocks where 
we extract useful hinds on the entropy generation leaving the study of more realistic shock 
for an upcoming paper. 

In this case the shock <3?(x + ,r) can be determined and the trapped surface ends at 
r = th with th determined from 

/iir F 
¥ , b\T H ) = — (4.5) 

with E ~ s5, and s the center of mass energy of the collision. The area of the trapped 
surface is 

r-r H 

Atrapped^ / b 3 (4.6) 
J oo 

We may therefore estimate the energy dependence of the trapped area for different bulk 
geometries. 

• For non-confining scaling theories b ~ r~ 7 , with 1 < 7 < 00. The AdS case corre- 

37-1 

sponds to 7 = 1. We obtain At ~ s 6 ^ . This agrees with previous estimates for the 
AdS case. 

• Confining backgrounds that are scale invariant in the IR, p6f| , with b(r) ~ (tq — r) s , 
5 > |. In this case we obtain At ~ s~es~ at high energy. In this case the exponent 
varies between I and 1. 

• Confining backgrounds with b(r) ~ e~( Ar ) a . In this case At ~ s2(logs) - ^ - at high 
energy. 

• Confining backgrounds with b(r) ~ e ( r - r o) . In this case A t ~ 52 (log s)~^~ at high 
energy. 

Note that in all cases the entropy production is larger than AdS. 

This simple context does not include the realistic structure of the transverse space 
distributions of real world heavy-ion collisions, but captures basic tenets of the approach. In 
particular it shows that the assumption of the gravitational description up to the boundary 
suggests a stronger s dependence of the multiplicity than seen experimentally. This is a 
consequence of string interactions of gravity up to the boundary. In a future publication we 
will analyse realistic collisions with localised energy distributions in the transverse space. 



5. Outlook 



The gravitational description of heavy-ion collisions, is promising to revolutionize the de- 
scription of both strong coupling physics in QCD and that of black hole formations and 
evaporation. We now have concrete and detailed experimental measurements of black-hole 
formation and evaporation that can help analyse in more detail the physics of (QCD) black 
holes, and give the only tool to calculate thermalization, hydrodynamics and black-hole 
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decay in heavy ion collisions. More analytical and numerical work is needed in order to 
analyze in detail these effects, but the technical tools seem to currently available. 

Already CMS has seen hints of collective behavior in high-multiplicity pp collisions at 
LHC. With higher energy, high-multiplicity pp data, we will have a new arena for 5 black- 
hole formation and evaporation. It is debatable whether the black holes in that case would 
be of the large kind. If they are "small" they will be even closer (but not identical) to the 
Schwartzschild -like black holes people expected from higher dimensional theories. 
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